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Abstract — In this paper we describe by a number of examples 
how to deduce one single characterizing higher order differential 
equation for output quantities of an analog circuit. 

In the linear case, we apply basic "symbolic" methods from 
linear algebra to the system of differential equations which is 
used to model the analog circuit. For nonlinear circuits and their 
corresponding nonlinear differential equations, we show how to 
employ computer algebra tools implemented in Maple, which are 
based on differential algebra. 
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I. Introduction 

Usually the input-output response of a linear time-invariant 
circuit is described in the frequency domain by its transfer 
function, i.e. a single rational function. This translates directly 
into a linear differential equation with constant coefficients in 
the time domain. The advantage of this approach is, that in any 
guise, only one (differential) equation is needed to completely 
describe the quantity a designer is interested in. 

This method fails miserably, when a transformation from 
the time to the frequency domain is not possible, e.g., when 
nonlinear circuits have to be examined, which lead to systems 
of nonlinear differential equations. Even though most of the 
times, these can be given in symbolic terms, any single 
quantity in the circuit usually is described by a "waveform", 
which results from assigning a numerical value to each 
symbol, followed by a numerical calculation using computer 
simulation. The advantage of having only one describing 
equation seems to have been irrevocably lost, when nonlinear 
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circuits are considered. Thus, up until now, nonlinear circuits 
seemed nearly inaccessible to most symbolic approaches. 

This problem is not a new one, and it is not limited 
to the area of analog circuits alone. Several years ago, re- 
searchers in nonlinear control theory have proposed to use 
constructive methods from differential algebra to tackle their 
problems. At the same time - and inspired in part by this 
proposal - mathematicians started to implement algorithms 
from differential algebra, which had akeady been formulated 
in the 1950s. These programmes became part of the MAPLE 
computer algebra system. 

At the SMACD-meeting in 1998 G. Carra FerrcQ gave 
examples of how to transform systems of nonlinear differ- 
ential equations containing certain transcendental functions, 
that arise from analog circuits, into systems of nonlinear 
"algebraic" differential equations uJ], and thus brought the 
area of constructive differential algebra in contact with the 
area of symbolic circuit analysis and design. In this paper we 
will take this approach several steps further and show, how 
single equations for any quantity in a circuit can be derived 
from these systems. We will be able to give a new and easy 
algorithm for linear circuits, which works in the time domain, 
and uses only differentiation and Gaussian elimination. For 
nonlinear circuits we will resort to the algorithms already 
implemented in MAPLE. We will comment on how to apply 
them and produce several examples. 

II. Linear Circuits and linear Systems of 
Differential Equations 

To get a first flavour of things, we start our discussion 
with linear circuits and their corresponding systems of linear 
differential equations. But, instead of working in the frequency 
domain, using the Laplace-transform to deduce the transfer 
function for a sought-after quantity and thus its characterizing 
differential equation, we will remain in the time domain. 

A. Linear State Equations 

In this section we will present a new algorithm, that, 
starting from a set of state equations, only uses repeated 
differentiation and, finally, Gaussian elimination, to compute 
the single differential equation for any given quantity. We will 

'Note added in 2008: Between the time of this article's first publication 
and the publication of the electronic version, Giuseppa Carra Ferro passed 
away on March 22nd, 2007. Thus a dedication to her memory was added to 
this electronic version. 



restrict our presentation to just three state variables xi,X2,X3, 
but it is easy to extend the method, shown below, to any 
number of state variables. 

So, let us suppose, that a given linear circuit can be 
described by a set of linear state equations. We ask for 
one single differential equation describing, without loss of 
generality, the state variable xi. Again it is easy to handle 
other state variables or any output variable y, which is a linear 
combination of state variables and inputs, in an analogous 
manner 

Let the system be given by 

xi{t) = an • xi{t) + ai2 ■ X2it) + 013 • x^it) + ei{t), 
X2{t) = a2i ■ xi{t) + 022 • X2{t) + a23 ■ x^it) + e2{t), (1) 
Xsit) = 031 • Xi{t) + 032 • X2{t) + 033 • X3{t) + 63(1), 

where xi, 2:2, 2:3 denote the state variables and 61,62,63 
represent linear combinations of the inputs (and eventually 
their derivatives). 

Clearly, if a;i,2:2,a;3 satisfy ([B, their derivative^ will 
satisfy 

xi{t) = an • xi{t) + ai2 ■ X2{t) + 013 • X3(t) + ei{t), 

X2{t) = 021 ■ Xi(t) + 022 ■ X2{t) + 023 • X3,{t) + 62(1), (2) 
Xsit) = 031 ■ Xi{t) + 032 • X2{t) + 033 • X3{t) + 63(1). 

If n > 3 state variables are given, we have to repeat the above 
procedure n — 1 times. It is a well known observation (which 
will be shown as a byproduct of our algorithm), that n linear 
first order state equations lead to one nth order differential 
equation for any single quantity. Thus, in our example, we 
need one further equation for the third derivative i'l of xi, 
which we get by differentiating once again the first equation 
of©. 

xi{t) = an ■ xi{t) + ai2 ■ X2{t) + 013 • x:i{t) + ei(t). (3) 

Next, we write down all of the above equations into one 
system, where the vector of variables is given by all the 
derivatives of all state variables, that have been produced by 
the above procedure. 
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We point out the fact, that the variables should be arranged 
according to the ordering X2{t) > X3{t) > X2{t) > X3{t) > 
X2{t) > X3{t) > 'i'i(i) > xi{t) > xi{t) > xi{t). For our 
algorithm to work, it is essential, that the last entries of the 
vector are the derivatives of the variable, for which we want 
to deduce the differential equation, in decreasing order 

^In the sequel wherever necessai'y, we will assume that all expressions are 
differentiable over a suitable extension field of the real numbers. 



The final step is to use Gaussian elimination to convert the 
system into one in upper triangular form 
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(5) 
where ^(t) is a linear combination of the functions ei{t), 
ei{t), 6i(i), e2(i), 62(i), e3{t), and e3(i). 

Since we had (n—l)-n+l (obviously) linearly independent 
equations for n^ + 1 variables, the last row will give one 
equation for the last n + 1 variables, which according to 
our special ordering are all derivatives of xi. Thus we have 
deduced the differential equation for xi, which we were 
looking for This is the same equation, which we would 
have found, if we had worked in the frequency domain and 
had calculated the transfer function for the Laplace transform 
C{xi) of xi. ■ 

B. Linear SemiState Equations 

It might happen, that a linear time-invariant circuit does not 
possess a description by state equations. Nevertheless it may 
be describable by so-called semistate equations, i.e. equations 
of the form 

Ex(<) = Ax(<) + Bu(i), y(i) = Cx(t) + Du(t), (6) 

in which A, B, C, D and E are constant matrices, E being 
singular, x(t) denotes the semistate vector, u(i) the vector of 
inputs and y{t) the vector of outputs. 

The algorithm given above can be adapted to this situation, 
but we have to keep in mind, that even if n semistate equations 
are given, the differential equation for any quantity might be 
of degree strictly less than n. We will see this effect in the 
example given below, which recently appeared in \^. 

Example 1: (OTA-Circuit) 
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The OTA circuit shown above can be described by the semi- 
state system 
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where x^{t) and x^{t) denote the node to ground vohages 
of the respective nodes and 0:3 (i) denotes the current into node 
®- 

As shown in |2], the transfer function H,{s) — 



given by 
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H{s) = gm3 



Cs + gmi 
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(7) 



We will deduce the corresponding differential equation for 
x^ [t) in the time domain, using the above system of semistate 
equations and a slight modification of the algorithm for state 
equations. This algorithm is closer to the one, that will be used 
for nonlinear circuits. It is based on two main principles: 

1) The variables Xf^^ ^(2)^ ^^ ^^'^ their derivatives, are 
supposed to be ordered by 

X3 < X3 < X3 < ■ ■ ■ < x^ < x^ < x^ < ... 

< ^(1) < ^® < ^CD < ■ ■ ■ 

When a term is to be eliminated, we always choose the 
term of highest order 

2) Equations are only differentiated when "necessary". 
So, let us get into details; the system was given by 

C x^ - C x^ = -g„i4 2^0 + .g,„4 x^ + X3 (8) 

-CXf^ + Cx^ ^-gjnsXQ (9) 

= ^Q-ei (10) 

The term of highest order appearing in (O - ( fTOb is x^^rs. 
Equation ( flOb gives 

2^0 = ei (11) 

This can be used to eliminate x^ from ([8]) and (|9|l. Further- 
more, after differentiating (fTTT i. we get 



x^=ei. 



(12) 



Thus we are able to remove all instances of x^ in the first 
two equations. Consequently we are led to: 

Cx^ = Cei-.gm3ei (13) 

Cx^ = Cei + grniSi - grniX^ - X3 (14) 

Clearly, these two equations imply the equality 
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which gives 
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Putting this into ( fT4] i. we arrive at 



C X3 
Cei - 3^361 = • {gm3 + 5m4)ei - C ■ , (18) 

9m,A 9mA 



which finally results in 

C ■ ±3= C ■ gm3 ■ ei + gva3 ' 9mA ' Gi. 



(19) 



This is the time domain equivalent of the transfer function ^, 
which we wanted to deduce. ■ 

III. Constructive Differential Algebra and the 

dif f alg-PACKAGE IN maple 

This is not the time and the space to give even a cursory 
treatment of those parts of differential algebra, which are 
needed to understand the sometimes subtle generalization to 
the nonlinear case of the algorithms shown above. For our 
purposes it is enough to know, that, cum grano salis, all the 
mechanisms are already visible in the example of the OTA 
circuit. 

Fortunately there already exist several implementations of 
the necessary algorithms within the computer algebra system 
MAPLE. We will show by way of the above example the 
workings of one of these, the dif f alg-package, created 
by F Boulier yt], ||j] and improved by E. Hubert |,5J et al. 
A more detailed description (suitable for beginners with a 
mathematical background) can be found in the world wide 
web iHl, H. 

After starting a MAPLE-session, one first has to load the 
package diffalg: 

> with(diffalg); 
[Rosen feld_Groebner, belongs Jto^ deltaJeader, . . .] 

(The second line in slanted notation represents the output 
produced by diffalg.) 

After initialization, one has to enter the set of differential 
equations under consideration. This has to be done in form 
of so-called differential polynomials. These are polynomials 
in the unknown functions xi, . . . , Xm, their (time) derivatives 
xi"' := D^-^^Xi := jt^Xi, 1 < i < m, a & N, the excitations 
ei, . . . , Cfc and their derivatives, again. 

In our example, we get three differential polynomials, which 
read in MAPLE-notation: 

> p_l := C*diff (x_l(t),t)-C*diff (x_2(t),t)+gjn4* 
x_l(t) - g_m4 * x_2(t) - x_3(t); 

> p_2 := -C * diff (x_l(t),t) + C * dif f (x_2(t), t) + 
g_m3 * x_l(t); 

> p_3 := x_l(t) - e_l(t); 

Next we have to tell the programme, which symbols it 
has to treat as constants. This is done with the command 



field-extension, dif falg assumes, that we work over 
the rational numbers as ground field, where any further 
constants are considered as lying in a transcendental field 
extension of Q (i.e. we are allowed to divide by constants 
different from 0, and constants do not satisfy any algebraic 
relations). If we work with symbols, e.g. for capacitors, 
resistors etc., this poses no problem. If we work with real 
coefficients (e.g. floating point numbers or algebraic numbers 
like \/2) major problems may arise. In our case we define 

> K := f ield_extension(trainscendental_elements = 
[C,gjii3,gjn4]); 

K :— ground_field 

Finally we define a so-called differential ring, which is 
supposed to contain all the objects of interest (i.e. differential 
polynomials and constants), and in which we are allowed to 
do the following operations 

1) multiply a differential polynomials with a constant; 

2) add and multiply differential polynomials; 

3) differentiate differential polynomials (if a constant is 
differentiated, the result is 0). 

As we have seen above, it is very important, that we define 
an ordering on differential monomials. This is needed to 
control the elimination process. For this purpose, dif falg 
asks for a ranking of the time dependent variables, from which 
it produces the obvious "elimination ordering". The variable, 
for which we want to know the differential equation, should 
be the last before the excitations. 
> R :=:dif f erentialjring(ranking = [x_l, x_2, x_3, e_l], 

derivations — [t], field_of .constants = K, 

notation = dif f ); 

R := ODE.ring 

The command Rosenf eld_Groebner lies at the heart of 
dif falg. It produces minimal sets of differential polynomi- 
als generating the differential polynomials, we have entered. 

> GE := Rosenfeld_Groebner({p_l,p_2,p_3},R); 

GE := [regular] 
GE is a list and may contain several component^ These 
components contain the sought-after differential equations, 
which can be listed with the help of the rewrite_rules 
command; in our example this gives 

> rewrite_rules(GE[l]); 

[x.m = e_l(i), 

-x_3(t) + e.l{t)g.m3 + e_l(t)g_m4 



x.2{t) 



gjiuA 



d 



g.m3 e_l(i)g_m4 + C ^e_l(t) 
^x.3(0^ ^ ^^* 



Mark, that the last entry is the differential equation for xs, 
which we had deduced in example [T] 

'This is one of the subtleties of working with dif falg, the details of 
which we do not want to present here. 



IV. Nonlinear Circuits and their corresponding 

DIFFERENTIAL EQUATIONS 

Now, that we know, how the MAPLE package dif falg can 
be used to deduce the time equivalent of the transfer function 
from any system of linear differential equations (with constant 
coefficients), in this section we will apply dif falg in an 
analogous manner to several systems of nonlinear equations, 
which come from nonlinear circuits. 

Example 2: (Damped resonant circuit with a nonlinear 
inductivity) 

As our first nonlinear example we have chosen a well known 
circuit from the literature, which leads to the Duffing equation 
(cp. f^, example 11-1, or fot], chapter 1.3.2.). 



Vo(t) 




The resistor and the capacitor are assumed to be linear, i.e., 
they are described by 

ic{t)=Cvc{t) and vnit) = RtR{t). (20) 

The inductor is assumed to be nonlinear, being described by 

VLit) = *(i), (21) 

where the current i^ (t) is approximated by the cubic 

tL{t)=a-il;{t) + b-i^{tf. (22) 

Finally, Kirchhoff ' s equation lead to 

init) = iL{t) = ic{t) and VR{t) + VL{t) + vc{t) = ~v^{t). 

(23) 
Equations (|20]|-(|23T| are translated into their corresponding dif- 
ferential polynomials and are used as input for the dif f alg- 
routine. This produces as part of the output of the subroutine 
rewrite _rules; 



dt 



sV-Ct) 



C 
CRa f^V'W) +iCRbi^{tf f-^^m 

C ' 

(24) 

which is formula 1 .30 in fv] and is an equivalent of the Duffing 
equation: 



Example 3: (Chua's Circuit) Our next example is Chua's 
circuit, the first example of a physically realizable circuit. 



showing chaotic behavior OlOll . This circuit consists of two 
hnear capacitors, one linear inductor, one linear resistor and 
one nonlinear resistor, as shown below. 

Vri 




Since it is easy to write down the equations for the linear 
elements and the Kirchhoff equations, we concentrate on 
the mathematical description of the nonlinear resistor The 
classic description of Rni, using piecewise linear functions, 
is unsuitable for our purposes, because it does not satisfy 
the right differentiability conditions. Thus we use the one, 
presented e.g. in 1.11.] . where the negative arctangent is used 
to produce the nonlinear v-i characteristic of i?„; . This leads 

Vnl{t) 



inl{t) 



-Iq ■ arctan 



Vn 



(26) 



We are now confronted with a new problem: the arctangent 
is a transcendental function, thus ( l26l l would not give a 
differential polynomial, as needed. One procedural solution to 
this problem was presented before in 1 1]. In the present paper. 



we are satisfied with the fact, that 
algebraic differential equation 



implies yet another 
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which obviously is given by a differential polynomial and 
which we can use as input for diffalg instead of ( |26] l. 
From this equation together with the Kirchhoff equations and 
the characterizing equations of the linear elements, diffalg 
produces the differential equation 

' (Ci + C2)L X + CiRx + X 



a;(4) = - 



o 1 C2 RL 



-lo- 
+ 2/0 

-In- 



Vo 



x^ + Vo^ 
Vo 



i^' + Vo'f 
Vo jGx^ - 2Vo^) 
(x^ + Vif 



[C2LRX +Lx + Rx] 

L XX ■ [3 C2R X + x\ 



C2LR'. 



(28) 



where x{t) denotes the voltage vci (t) through the capacitor 
Ci- It has to be said that the final result given by diffalg 
looks slightly different, since it is given in expanded form, 
i.e. the numerator consists of 31 summands. Formula ( l28T l has 
been reached at, only after some laborous post-processing. ■ 



Example 4: (Simple Model of a Peak Rectifier Circuit) 

Now we are going to show how to handle diodes (and 
consequently by way of the Ebers-Moll model the large signal 



behaviour of BJTs) in electric circuits. As an example we have 
chosen a simple model of a peak rectifier circuit as seen in 
Iil2,1 , chapter 3.7 pp. 185ff. 
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Again we concentrate on the only nonlinear element in the 
circuit - the diode. It is well known, that the v-i characteristic 
of a nonideal diode can be approximately described by 

'voity 



ioit) 



exp 



Vt 



1 



(29) 



where Ig is the saturation current and Vr is the thermal voltage 
- quantities, which we consider constant during the course of 
our analysis. 

As before we have to translate a transcendental equation 
into a differential polynomial. This can be done easily by dif- 
ferentiating ( |29] l once, which due to the chain rule ^ ioit) 

dvn 



D ' li'^D (t) results in the equation 



d , , 1 f d 



[tD{t)+Is]. 



(30) 



This time diffalg produces the following second order 
differential equation for the output voltage Va{t) : 

.. _ Ra- [VrVa - (Vo-Va) ■ Vx] + [Va ' ReVx] 



CRa ■ {VrRa 

where we have set Vx ■= {CRaVa - 



-ReVx) 

Va + Rah) ■ 

For the purpose of comparing this result to that appearing 
in the literature, we give the differential equation in case of 
an ideal voltage source, i.e. R^ = 051. It is given by 

.. 1 / . 1 



[Vo - Va\ 



[CRaVa 



Rah 



V. Further Examples 

Example 5: (Diode Circuit with LC-Load and Nonideal 
Voltage source) 
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The above circuit may not be of much practical interest. 
Nevertheless it is a good test of the power of our approach 
(and the capabilities of diffalg), since it slightly generalizes 
example |4] and we increase the number of dynamic elements 
in our circuit. Again, the diode is assumed to be nonideal, 
given by (|30] |. 



!'„(<) 



The output voltage x{t) :— vc{t) is described by the 
equation 

r 1 . x + cLi 1 

R 



X — vo Vt ■ CL 



L 



[x + CLx) + {x — vq) 



Vt {x - Vq) 



(33) 



In the course of our investigations, we have tried a number 
of larger circuits, which in principle are accessible to our 
approach. We met two main obstacles, which are natural in 
the "business" of symbolic methods: 

1) the combinatorial explosion, resulting in a larger and 
larger number of terms contained in the final differential 
equation, and 

2) the massive increase in time, needed by diffalg to 
produce this equation. 

In the sequel we give a short report on these experiments: 
Example 6: (Peak Rectifier with Power Transformer) 




This kind of circuit is described in the introduction to 



chapter 3.7 in 111211 . As shown above, we have used the model 
given by Horneber in his PhD-thesis 11311 . section 14.1. There, 
it is the smaller of two examples, the other being the "Ring 
Modulator", which has become a benchmark in the numerical 
analysis of initial value problems il4tl . 

From our point of view, we can tell, that diffalg, 
although needing substantial more time than in the previous 
examples (1-2 minutes instead of only seconds), is able to 
produce a fifth order differential equation for the output volt- 
age Va- Unfortunately, we are not able to reproduce this result 
here, since the initial output even after some simplifications 
contains more than 600 summands. Thus some more "post 
processing" is needed to get an intelligible result. Even with 
the help of other facilities of MAPLE this work has not been 
finished, yeQ ■ 

Finally we have tried our approach on a "simple" single-stage 
common-emitter amplifier (" 1.121 . chapter 4.11) as modelled 
in lll5ll and on the above mentioned ring modulator of 
Horneber In both cases, up until now, even though we have 
used several days of computing time, we were not able to 
produce any results. Although the latter - very ambitious 
- example (which presumably will lead to an differential 
equation of order 18) may be beyond the scope of any 
computer algebra system for some time, the first should be 
within our grasp and should be attacked further. 

'*Note added in 2008: Meanwhile tliese calculations have been done. The 
end result still is to unwieldy to be presented here. Furthermore the collecting 
and combining of fully symboHc terms by hand has turned out to be so error- 
prone that, in the opinion of the author, some kind of additional "plausibility 
measures" need to be introduced. 



VI. Conclusion 

In this paper we have shown, how, using constructive meth- 
ods from differential algebra and one of their realizations - the 
package diffalg of the computer algebra system MAPLE - 
linear and nonlinear circuits can be described by a single 
differential equation. In the future it will be necessary to 
further examine the power of this approach, i.e. to find more 
and larger circuits, which can be treated this way. Furthermore, 
if the number of these circuits is large enough, methods have to 
be found, that allow a fast and "easy" analysis of the resulting 
equations, analogous to the analysis of linear circuits by way 
of their transfer functions. 
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Note added to the Electronic Version 

In this electronic document, some small typographical errors 
of the printed version were corrected. This especially refers to 
formulas ([T8]l and ( fT9] l. 

Furthermore, for the convenience of the reader the abstract 
has been rewritten, and keywords, an MSC classification, 
and a short CV according to IEEE standards have been 
added. URLs have been checked again, and, where neces- 
sary, have been updated. Finally the dedication has been 
expanded. The main body of the article, however, remains 
unchanged. (April 17th, 2008) 



